Observation of stable phase jump lines in convection of a twisted nematic 
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Abstract 

We report observations of stable, localized, line-like structures in the spatially periodic pattern formed 
by nematic electroconvection, along which the phase of the pattern jumps by tt. With increasing electric 
voltage, these lines form a grid-like structure that goes over into a structure indistinguishable from the well 
known grid pattern. We present theoretical arguments that suggest that the twisted cell geometry we are 
using is indirectly stabilizing the phase jump lines, and that the PJL lattice is caused by an interaction of 
phase jump hnes and a zig-zag instability of the surrounding pattern. 
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1. INTRODUCTION 



Usually, topological defects in the phase of periodic spatial modulations have co-dimension 
two. They are point-like in two dimensions, line-like in three dimensions, and in one dimension 
they reduce, for dynamical systems, to short events, localized in space and timeiJ]. Here we report 
the observation of stable, self-organized phase defects that extend to line-like objects in 2d: Phase 
jump lines (PJL) stretched parallel to the wavevector of the modulation, characterized by a change 
of the phase of the spatial modulation by tt on a short distance (Fig[T]). 

The phenomenon was observed in the periodic pattern of convection rolls forming in thin layers 
of the nematic liquid crystal MBBA (4-methoxybenzylidene-4'-butylaniline) under the influence 
of an electric ac field (p) that is oriented normal to the {x-y) layer and the alignment of the nematic 
director. Electroconvection (EC) in nematic liquid crystals has been studied extensively. After 



early studies by E. Dubois- Violette and others ||2|] aimed at understanding the basic mechanism 
of convection, EC later became popular as a model for pattern formation in anisotropic systems. 
Depending on the experimental conditions, the convection rolls are aligned either preferentially 
normal to the symmetry axis (x) of the system — defined by the direction at which the nematic 
director is anchored at the boundaries — or oblique at some preferred angel. In the oblique-roll 
case, two convective modes related to each other by a reflection at the x axis coexist jsi ^. Near 
the onset of convection, descriptions in terms of one (normal rolls) or two (oblique rolls) time- 
dependent Ginzburg-Landau equations could be established [5]. These equations are of potential 
type and predict a relaxation to a simple equilibrium state. 

The more puzzling was the rich dynamics of the system observed for stronger electric fields, 
with scenarios such as those named "fluctuating Williams domains", "chevron patterns", "abnor- 
mal roUs", "grid patterns", or "defect lattices" being reported lalllalalialiillillillilliy, 
[l^lisl- Only in recent years it became clear that most of these observations can be understood 
by the non-relaxational interaction of the convective mode(s) with only one additional, weakly 
damped mode that describes an in-plane rotation of the nematic director. A Ginzburg-Landau- 
type model for the non-relaxational coupling of these two (or three) modes, called standard model 
below, explains most of the experimentaly observed scenarios at least qualitatively, and predicts a 
few more Ofl. 

There are several indications in the literature that hint at the existence of stable PJL in nematic 
electroconvection. For example, for the defect-chaotic regime of normal rolls called "fluctuating 
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Williams domains" it is well known that the "shape" of topological defects, particularly near the 
events of defect creations and annihilations, can be highly anisotropic. Their extension in x di- 
rection is considerably larger than in y direction. As such this deformation could be explained 
by the anisotropy of the linear dispersion relation of the convective mode, which can be removed 
by a simple rescaling of the x axis. But in some cases the defects stretch considerably and are 



then more naturally characterized as "phase jump lines" ||21 



2311 . These structures seem to 



be related to weakly unstable saddles corresponding to the stable PJL reported here. Stable PJL 
have also been seen in numerical simulations of the standard model j^l- Yet, a general theory is 
missing, and it remains unclear if or under which conditions stable PJL would be observable in 
electroconvection . 

In our experiments, the conventional EC setup was modified by twisting the nematic within 
the layer: Instead of anchoring the nematic director at the surfaces of the layer parallel to the x 
axis, the anchoring direction is rotated by a fixed angle in opposite directions at the two surfaces 
iHQHQ- This geometry seems to be favorable for PJL. 

In the following we first describe the details of the twisted EC setup, and describe the observed 
PJL and structures they form. We then discuss some theoretical ideas that understanding these 
observations. 



2. EXPERIMENTAL SETUP 

The nematic liquid crystal MBBA doped with 0.1wt% tetra-n-butylammonium bromide to sta- 
bilize the conductivity was sandwiched between two glass plates covered with transparent indium 
tin oxide thin-film electrodes. The lateral extensions of the liquid crystal layer were 2cm x 2cm 
and its thickness d = 50/im. The surfaces of the plates were coated with polyvinyl alcohol and 
rubbed in order to attain a planar anchoring of the nematic molecules. Temperature of the cell was 
controlled at 25±0.01 °C. A standard wave-generator was used to apply ac voltages V = C(IOOV) 
of frequencies / up to several kHz to the electrodes. The structure of the cell is shown schemat- 
ically in Fig. [T] The anchoring direction for each plate is represented by an arrow. By rotating 
the anchoring direction by an angle a = 7r/4 out of the x axis at the lower plate, and in opposite 
direction by the same amount at the upper plate, the nematic director is twisted along the z-axis in 
the ground state [40] (Fig.[T]). 
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3. EXPERIMENTAL RESULTS 



The control parameters of the system are the externally applied voltage and its frequency. A 
phase diagram of the observed patterns is shown in the Fig. |21 The cutoff frequency where 
the conductive mode of electroconvection is superseded by the dielectric mode, was above 1 kHz. 
Measurements were performed for / < fc, i.e., in the conductive regime. As the voltage was 
increased, normal rolls with a wavevector parallel to the a;-axis (at an angle of a to each anchoring 
direction) were observed. Besides the well-known states including normal rolls (parallel rolls), 
grid pattern (rectangular cells), and dynamic scattering modes (3D turbulent states), we found 
stably aligned phase jump lines (PJL) in a wide frequency range of the conductive regime. 

Figures |3^-c show the shadow-graph image of PJL observed under a polarization microscope 
(polars II x). PJL can have finite length. At the endpoints of PJL the phase difference between the 
adjacent domains gradually decreases to zero. As the voltage is increased, the typical distance L 
between PJL in y direction decreases. 

In order to determine L, which is below expressed in terms of the wave number K := 27i/L, 
quantitatively, two methods were employed. The first makes direct use of the Fourier transforms 
of recorded images (Figs.|3j5',c'). The other method is to count the number n ll4li of PJL along 
a line of length E parallel to y, which gives K = 27m/ E. For high voltage, the Fourier analysis 
was preferred. At lower voltages we used the counting method, since the small number of PJL 
would not lead to sharp peaks in the Fourier transform. We verified that in the intermediate range 
the results of the two methods coincide. 

In Fig. m the measured wave number K is plotted against the normalized distance e = (K^ — 
V^)/V^ from the threshold of convection Vc. The data indicate that K goes to zero following a 
square-root law 

K = c{e-e^Y/' (1) 

with ei > 0. For / = 400Hz we obtain c = 50.5mm~^ and ci = 0.22 and for / = 600Hz data is 
well described by c = 95.3mm^^ and ei = 0.3L 

Some observations indicate that there are at least two different regimes of PJL. At low e, the 
spacing between PJL is very irregular (Fig. |l^,b) and the PJL continue to move, although very 
slowly. At high e, where the spacing between PJL is smaller, it becomes also more regular (as 
can be seen from the sharp peaks in the Fourier transform of Fig. |3t shown in Fig. |3t'), and a 
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steady state is reached. We refer to the resulting structure as a PJL grid. The onset voltage of 
PJL grids is well above the onset of PJL (Fig. |2l). At this point the K vs e relation also starts 
to deviate significantly from the square-root law valid for smaller e, as shown in Fig. HI The 
deviation is towards larger K, i.e. towards denser PJL, contrary to what one would expect from 
a repulsive interaction between PJL. The PJL grids are visually indistinguishable from the well- 



known grid pattem ll29l] consisting of two overlapping sets of oblique rolls. If the two structures are 
indeed identical this means that we found a new route from normal rolls to oblique rolls, different 
from the conventional transition via a zig-zag instability. But further investigations are required to 
confirm this hypothesis. For the low-e range, on the other hand, some first steps towards a better 
understanding have been made. 

4. THEORY 

A clear theoretical picture of PJL, including the conditions for their existence as well a descrip- 



tion of their interaction with the surrounding roll pattern and with eac 



1 other, has not yet emerged. 



Even though other mechanisms leading to stable PJL appear to exist 113011 . the simulations of Zhao 
strongly suggest that the observed PJL can be described by the standard model for normal- 
roll electroconvection. We shall here investigate three aspects related to this problem that help 
understanding why PJL stabilize in twisted cells and how they are forming PJL grids. 
The standard model is given by 

TdtA =i(3AdyCt> + (e + Cdl + 

^ . (2a) 



^^dt<P = -H + K^dl<P + K^dl(t) 

(2b) 

+ 2T\m{A*{dy-iqc<p)A}. 

It describes the interaction of two spatio-temporally slow variables: The complex amplitude A of 
the convection pattern ~ (Ae^''"^ + c.c), and the amount (p of in-plane rotation of the nematic 
director. By the requirements of structural stability and according to theoretical and experimental 
results for the untwisted geometry, the time-constants r and 71, the coherence lengths ^y, the 
critical wave number q^, the geometry factors Ci and C2, the Landau coefficient g, the damping 
constant h, and the two elastic constants Ki and K3 are all positive. The sign of P depends on the 
applied frequency and the experimental geometry and is difficult to determine a priori. The value 



of r is always found to be negative and we shall assume this hereafter. It is easily seen that with 
negative T convection rolls {A ^ 0) have a destabilizing effect on 0. This effect is crucial for most 
of the complex patterns observed in electroconvection, and PJL are no exception. 



4.1. Stabilization of PJL 

After scaling out the natural e dependence of A, and length and time scales (d^ dy A r~~> 
~ e^/^, dt ~ e), the standard model can be reformulated by the following, equivalent system 
of PDE: 

TdtA =(1 + 92 + dl + 2i C^(Pdy - C202 - l^n^ + il3Ady(P, (3a) 



y 

I 



d^(t) = -h(t) + dl(l) + K^dlct) - 2\T\ \m{A*{dy - icfyA}. (3b) 



The coefficients and variables in system dSj) are not identical to those denoted with the same sym- 
bols in system ©, but are related to these by simple scaling transformations. For all coefficients 
but h this relation is independent of e. The coefficient h in Eq. (l3bb scales as e^^, i.e., h plays the 
role of the main control parameter in the rescaled equations ©. As /i decreases, the first instability 
of the homogeneous solution A= 1, = is to abnormal rolls at/i = 2|r|if/3<0 and to zig-zag 
modulations at /i = 2|r|(l + /3) if /3 > 0. 
The system also has a steady-state solution 



0, A = tanh (y/Vt^ (4) 



that describes PJL. For large h, i.e. close to the threshold of convection, PJL are unstable. This can 
be seen by considering the limit h ^ 00 where = is fixed. Then solution Q is unstable to the 
linear mode 6A = i sech(y/ ^/2) with growth rate I34I1 . However, by lowering h the additional 
degree of freedom can stabilize the PJL by an interaction through the last term in Eq. (l3bb . For a 
discussion of this process, we linearize the system Q with respect to perturbations of solution Q 
of the form 5 A = [ar{y) + iai{y)] exp(crt), 50 = f{y) exp(crt), which leads to 

-(jr + 1 -Stanh^^ + 9j 1 a^, (5) 



V2 
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and 



= (^-V2Ci sech^-^ + P tanh dy^ f 



y 

—a T + sech^—= + dy ] a. 
v2 



(6a) 
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-a-h- 2|r|tanh2^ + 5j j / 
( \/2 irlsech^^ - 2 irl tanh -^dy 



(6b) 



We are looking for eigenfunctions of these equations which are bounded for y — > ±oo. Equa- 
tion ^ is decoupled from Eqs. (I6al6bl ). Apart from the translational mode = dy tanh(y/-\/2), 
Eq. (jSj) has bounded solutions only with ra < —3/2 and does not contribute to the linear stability 
problem. The y axis can be divided into an inner region ||/ 1 < 5 near the PJL, and two outer regions 
where the system Q simplifies. For y > 0, for example, the substitutions sech — 0, tanh 1 
lead to 

0=(-ar + d^y)ai + (3dyf, (7a) 
0=(-a + 2|r|/i + 9j)/-2|r|9,ai, (7b) 

where we introduced the abbreviation fi = 1 — h/2\T\. In general, the eigenmodes of Eq. Q 
have to be calculated numerically. For example this can be done by "shooting" solutions start- 
ing at y = with aj = 1, f = fo, a'- = f = into the y > outer region such as to match 
ai{yi), a-(?/i), f{yi), f'{yi) for sufficiently large yi with a bounded outer solution. We shall here 
only be interested in the cases where the PJL stabilize prior to any destabilization of the homo- 
geneous state, i.e., for the cases that /i < — /3 and ji < 0. The bounded outer solutions are then 
generally given by linear combinations of the two solutions of dTJ of the form a^, / ~ exp{\y) 
with A < and 

A^ =^ + 0{a') (8) 
p + H 

or 

A2 = _2|r|(/3 + /i) + 0(a), (9) 

respectively. The value of ji where the PJL stabilize can be found by considering the limit cr — 0+ 
of the systems ^ and Two particularities of this limit will be used here: (i) The outer solution 
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corresponding to (EJ) goes over into a simple phase shift (/ = 0, = const.). Thus, / contains 
only contributions from the component corresponding to ^ and must satisfy 



dyf = -^-2\T\i/3 + fi) f. (10) 



(ii) By Eq. dTal) the quantity 



u = dya, + (3f (=0) (11) 



is conserved along y in the outer region [24]. When fi < —(3 and (flOb holds, this implies m = 
for bounded solutions. Equations (IIOII II) are necessary and sufficient conditions for the inner 
solutions to be bounded when continued io y oo. When the shooting method is used, the two 
conditions can be satisfied by adjusting the initial value /o of / at y = and the value of the 
control parameter /i. 

With Ci > 1 we could find neutrally stable modes of PJL with fi < —[3 for all parameter sets 
tested, in particular for positive as well as negative (3. For example. Fig. © shows the neutrally 
stable mode of the PJL for parameters V = —0.6, P = 0.3, Ci = 1.1, and r = 0.5. These 
parameters are suggested by numerical calculations for similar systems IbiIi . For this case we 
verified that there is no unstable mode, i.e.. Fig. Q shows the critical mode of a PJL. (For smaller 
|r|, e.g. with r = —0.06, there is another mode with positive growth rate. Then the feed-back 
loop through is too weak to stabilized the PJL.) 

As Ci is lowered, the critical /i approaches — /5 and the decay of the / component in the far 
field becomes flatter, as described by Eq. (flOb . Thus increases/decreases approximately linearly 
over a wide range according to (fTTt . For /i — > — /3, one expects the critical mode to converge to an 
unbounded secular solution with constant dyai and / in the outer range. 

It turns out that for Ci = 1 the linearized equations © have the analytic solution 

/ = 1, = (1 — /i) + /iy tanh (12) 

v2 

For /i = — /3 this solution satisfies conditions (flOb and dTTT) and is of the expected form. Combined 
with our numerical calculations, we conclude that the critical /i reaches — /3 exactly when Ci = 1. 
Thus, for Py > 0, the PJL can stabilize prior to the zig-zag instability at /x = — /? only if Ci > 1. 

For P < the first instability of the homogeneous state and the far field in system ^ is towards 
abnormal rolls at yU = 0. In contrast to the zig-zag instability at /i = — /3, the outer equations © 
have for = no secular solutions with a = 0. Instead, the solution space corresponding 
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to A = [Eq. (El)] is spanned by two linearly independent constant solutions, the pure phase 
shift {ai, f) = (1,0) and the pure director rotation (ai, f) = (0,1). Since both are bounded, 
bounded solutions of © for fi = are obtained whenever the contribution to the outer solutions 
corresponding to Eq. ® with A > vanishes. When the shooting method is used, this can be 
achieved by adjusting a single parameter, e.g. /q. Thus, at the abnormal-roll instability fi = 0, a 
neutrally stable linear mode of the PJL can always be found. 

By the interaction of / and in the core region, a certain linear combination of the constant 
outer modes is singled out, characterized by a specific ratio r = a^/ f (for y +oo). When 
Ci = 1, for example, the analytic solution (fT2ll gives r = -\/2. A perturbative analysis of the outer 
equations ® for small cr ~ 9^, ~ ;U under the constraint a,// = r + C(/i) yields a solution with 
Oj, / ~ exp(As?/) (as y +oo) where 



A. = , ,J +0(/i^) (13a) 



2rr|r 




2r2 


r 


T 





with the growth rate 

When As < this solution is bounded and an eigenmode of the PJL. Since a/\s = (3 /rr + 0{ij?), 
this mode can be stable only for r < 0, which is the case only for values of Ci much smaller or 
much larger than one. Usually one would assume r > 0. With — 2|r| r < (3 < 0, for example, 
([T3t then describe an unstable linear mode that becomes neutrally stable as /i approaches zero 
from below. 

In conclusion, these considerations show that within the framework of the model Q a stabi- 
lization of PJL of the form Q prior to any destabilization of the homogeneous state — as we seem 
to observe it in the experiments — is possible, but likely only for Ci > 1 and /3 > 0. 



4.2. The effect of the twisted anchoring on Ci 



In order to test if the twisted anchoring used in our experiment can cause a higher value of Ci, 
we analyzed a model for the z dependence of director twist 0(z) and convection mode A{z) and 



their interaction ||3C , 



35h . Keeping only terms that are relevant for determining Ci, the model reads 
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in appropriate dimensionless units 



TdtA 



A (14a) 



u-{P-4>Y + di 

dt^=dl^ + 2T\A\^{P-^). (14b) 

The coefficient u in Eq. (I14al) is a control parameter that depends on the applied voltage. P 
corresponds to the angle of the wavevector of a plane- wave convection mode to the x axis and the 
term — (P — 0)^ in Eq. ( I14al ) describes the preference of the convective modes to align normal 
to the director. The last term {d^) describes the diffusion of hydrodynamic variables along z. In 
Eq. (I14bl) the first term describes the stiffness of the in-plane director to twist and the second 
term the torque generated by misaligned rolls on the director. We assume layers of thickness 2 
without loss of generality and impose the boundary conditions ^(±1) = and, corresponding the 
twisted anchoring, 0(+l) = h, 0(— 1) = —h. The isotropy of the system is broken only by the 
boundary condition. In the bulk, isotropy is preserved, which is reflected by the equivariance under 
P^P + 5,(j)^(j) + 5. The model can be derived analytically in the limit that the wavelength 
of the roll pattern is much smaller than the layer thickness IboIi . However, since this condition is 
not satisfied for our experimental system, we expect the model to describe the experiment only 
qualitatively. 

The coefficient Ci entering the 2d description @ can be obtained by the following method. 
Noting that in the basic state {A = 0) the director relaxes to = bz, and that the fundamental 
mode of linear excitations of is 0i = cos(7r2;/2), we assume a weakly excited (j) = bz + K(j)i 
and search numerically for the threshold value Merit of the control parameter u at which the lowest 
eigenmode A{z) = Acnt{z) of Eq. i\4ah becomes critical {dt = 0). The orientation P of the 
roll wavevector is chosen such as to minimize Merit- This corresponds to letting A = exp(iPy) 
in Eq. (l3al) with P = Cicf). Following general prescriptions (e.g. \y\ ). the value of Ci is 
then given by the ratio of projections of two components of the r.h.s. of Eq. (I14bl) onto its adjoint 
eigenvector (= 0i) 



1 

2; 



^i2T\A\^Pdz 

Ci = (15) 

01 2T\A\'^4)dz 



1 



in the limit K 0. 



Without twist (b = 0) we reproduce the analytic result ||42|] Ci = 256/(277r2) = 0.96 



With increasing twist b, the value of Ci increases, until it reaches a maximum Ci = 1.14 at 6 = 8.1 
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(Fig. |6l). The result is independent of all other system parameters. Due to the purely geometric 
nature of this effect, it is plausible to assume that the twist leads to an increased value of Ci also 
in our experimental system. 



4.3. Coupling to zig-zag modulations 

Simulations 1241 show that PJL can also coexist with zig-zag modulations of the roll pattern. 
The PJL are then located at preferred phase angles of the zig-zag modulation. Thus, the distance L 
between neighbouring PJL will be controlled by the wavelength A of the zig-zag modulation. Even 
though they are not clearly visible in the experiments, small zig-zag modulations are probably also 
controlling the distance between PJL in PJL grids. 

To test this hypothesis, the wavelength of zig-zags is estimated. For an easy comparison with 
experiments, dimensional units are used in this calculation. In the absences of strong variations of 
the pattern amplitude \A\, the dynamics of the roll pattern can be described in terms of the phase 
9 = arg A (mod 27r) of the pattern alone. Together with the coupling to (p, this is description 
becomes tt32 



rdtO =ildle + + (3dy(t) + r(t)''dy<j), (16a) 
7i5t0 =2|r|(e - eAR)0 - g'f + K:,dl<P + K^dl^p - 2\V\edye (16b) 

with positive r, g, the threshold of abnormal rolls eAR = h/{2\V\) and all other parameters identical 
to those entering Eq. Q. While the derivation of the nonlinear terms in this system hinges on the 
assumption of that spatial and temporal variations of and are slow on the scales ix,y^^^ and 
re^^, which is strictly satisfied only for small e — eAR and (3, the linear part is valid uniformly in 
e < and (3, as long as the Ginzburg-Landau description Q is valid. 

For stationary solutions that depend only on y (dt, dx = 0) Eq. (I16al) can be integrated to 



ildyO + 13^ + r0V3 = J = const. (17) 

In the simple case that rolls are not tilted on the average {{dyO)^ = 0) the integration constant 
J vanishes. Using Eq. (fTTl to eliminate 9 from Eq. (I16bl) then yields an equation for of the 
form [|39|] 

K,d'y4> = -2|r|((l + /3/ey)e - eAR)0 + g\e)4>\ (18) 
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While the wavelength of zig-zag modulations can become arbitrary large when the nonlinearity in 
Eq. (fTSt becomes effective, the wavelength of small- amplitude modulations is sharply determined 
as 



A =27r, 



2|r|((l + /5/e2)e-6AR) 

— (19) 

Kiezz 



2 T eARie - ezz) 



where ezz = ^ar/ i'^+P/^y) denotes the threshold of the zig-zag instability. In order to estimate the 
value of 2|r|eAR/-ft'i, recall that the linear part of Eq. (I16bl) holds also for e = (no convection). 
In this case the relaxation of small in-plane rotations of the director is easily described by the 
linearized Leslie-Erickson equations. For the untwisted case, for example, they reduce to 

lidtuy = ihsdl + kudl + k22dl)ny. (20) 

Considering only the slowest mode of the system for fixed boundary conditions dX z = ±(i/2, 
Uy = (f) cos{7r z / d) , and projecting onto the corresponding adjoint cos{7r z/d), yields 

^A<l> = ihsdl + knd'y - k22^)(t>- (21) 

Comparing of Eqs. (fT6bll and (EB shows that 2|r|eAR/i^i = n'^k22/d'^kn. For MBBA at 25°C, 
for example, the ratio kii/k22 is ~ 1.6. For the twisted geometry one would expect some mixture 
of kii, k22, and ^33 to enter Ki instead of just kn, but since in the projection onto the adjoined 
mode the central part z ^ dominates, and there the director is still aligned parallel to x, these 
corrections are presumably small. Assuming a relation L = mA between zig-zag wavelength A 
and PJL distance L, one obtains the result 



L =2dm\ 



kiiezz 



k22ie-ezz)' ^^2) 



^ K -— / ^22 e - tzz 
dm V kii ezz 

Even when ezz is fitted from the measured data for K, Eq. (E^ makes a nontrivial prediction 
because ezz appears at two different positions in the formula. 

Using the fitted experimental curves, we find that Eq. (l22|l holds with m = 2.1 for the 400 Hz 
data and with m = 0.94 for the 600 Hz data. That is, in both cases we find a value of m of order 
one. In view of the coarseness of the theoretical estimates, and allowing for some experimental 

12 



error, we conclude that the observations fully supports the hypothesis that the lattice constant of 
the PJL-grid is controlled by small zig-zag modulations. The mechanism by which the amplitude 
of zig-zag modulations or, equivalently, the zig-zag wavelength is kept small remains unclear, 
however. 

5. CONCLUSION 

In this study, we reported the observation of stable PJL and described a process by which these 
PJL from a lattice that finally evolves to a structure indistinguishable from the well know grid 
pattern as the applied voltage increases. We argued that the stabilization of PJL is presumably an 
indirect effect of the twisted geometry of the nematic in the basic state: The twist leads to values 
of the coefficient Ci > 1, which in turn enables the stabilization of PJL prior to a destabilization 
of the roll pattern. 
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FIG. 1: Top: rubbing direction and the anchoring direction of nematic. Bottom: a typical phase jump Une 
at 400Hz, 13.10 Volts. The white bar corresponds to 100/xm. 
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FIG. 2: Phase Diagram of our experiments. Solid circles represent the onset voltage of convection rolls 
(Williams Domain); open square the onset voltage of PJL; open circle represents the onset of the PJL grid 
pattern (GP); solid square the onset voltage of the dynamic scattering mode (DSMl). The inset shows the 
part of the phase diagram that we discuss in this paper. The points a, b and c in the inset correspond to 
Figs.|3^-c. 
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FIG. 3: The upper row shows images of convection patterns with PJL, the lower row the modulus of the 
corresponding 2D Fourier transforms. These three pair of figures represent the transition from PJL to GP. 
Fig. a corresponds to 13.93 V (e = 0.47), Fig. b to 15.91 V (e = 0.91), and Fig. c to 19.23 V (e = 1.79), all 
at 400 Hz. The length of the scale bars is 200/im. 
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FIG. 6: The geometry factor Ci in dependence of the twist b for the simple model (IT4li . The dash-dotted 
line at Ci = 1 is a guide to the eye. 
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